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Joule heating in nanowires 
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We study the effect of Joule heating from electric currents flowing through ferromagnetic nanowires 
on the temperature of the nanowires and on the temperature of the substrate on which the nanowires 
are grown. The spatial current density distribution, the associated heat generation, and diffusion 
of heat is simulated within the nanowire and the substrate. We study several different nanowire 
and constriction geometries as well as different substrates: (thin) silicon nitride membranes, (thick) 
silicon wafers, and (thick) diamond wafers. The spatially resolved increase in temperature as a 
function of time is computed. 

For effectively three-dimensional substrates (where the substrate thickness greatly exceeds the 
nanowire length), we identify three different regimes of heat propagation through the substrate: 
regime (i), where the nanowire temperature increases approximately logarithmically as a function 
of time. In this regime, the nanowire temperature is well-described analytically by You et al. 
[APL89, 222513 (2006)]. We provide an analytical expression for the time tc that marks the upper 
applicability limit of the You model. After tc, the heat flow enters regime (ii), where the nanowire 
temperature stays constant while a hemispherical heat front carries the heat away from the wire and 
into the substrate. As the heat front reaches the boundary of the substrate, regime (iii) is entered 
where the nanowire and substrate temperature start to increase rapidly. 

For effectively two-dimensional substrates (where the nanowire length greatly exceeds the sub- 
strate thickness), there is only one regime in which the temperature increases logarithmically with 
time for large times, before the heat front reaches the substrate boundary. We provide an analytical 
expression, valid for all pulse durations, that allows one to accurately compute this temperature 
increase in the nanowire on thin substrates. 



I. INTRODUCTION 

Recently, there has been much interest both in fun- 
damental studies of spin torque transfei'^^ and in ef- 
forts to realize devices such as th e ra ce-track memory 
exploiting the spin torque transfer.!^ In either 
the present very large current densities have to be used 
to move domain walls and, more generally, to modify 
the ferromagnetic patterns. Associated with these large 
current densities is a substantial amount of Joule heat- 
ing that increases the temperature of the sample or de- 
vice. It is a crucial question to understand how strongly 
the temperature increases as this can affect the observed 
physics considerably, e.g. Refs. IMTUl and may even lead 
to a temporary breakdown of ferromagnetism if the Curie 
temperature is exceeded. The depinning of a domain wall 
could be due to a strong spin-current torque transfer, or 
as a result of the extreme heating of the material due to 
reduced magnetic pinning at elevated temperatures, or 
due to the intermittent suppression of ferromagnetism. 

You et alW^ have derived an analytical expression that 
allows to compute the increase of temperature for a 
nanowire (extending to plus and minus infinity in y- 
direction) with height h (in z-dircction) and width w (in 
x-direction) . The nanowire is attached to a semi- infinite 
substrate (which fills all space for z < 0) . The heating is 
due to an explicit term S{x, t) which can vary across the 
width of the wire and as a function of time. 



Meier et al^ use energy considerations to estimate the 
total amount of energy deposited into the nanowire and 
substrate system to show that - for their particular pa- 
rameters - the heating and associated temperature in- 
crease stays below the Curie temperature. 

In this work, we use a numerical multi-physics simu- 
lation approach which allows to determine the temper- 
ature distribution T(r, t) for all times t and positions 
r. Starting from a given geometry and an applied volt- 
age (or current), we compute the resulting current den- 
sity, the associated heat generation, and the temperature 
distribution. While such a numerical approach provides 
less insight than an analytical approximation, it allows 
us to exactly determine the temperature distribution for 
nanowires of finite length, nanowires with constrictions 
and thin substrates for which the assumption of an infi- 
nite thickness is inappropriate. Thus it evinces the limits 
of applicability of analytical approximations. 

While the simulation and analytical techniques used 
and developed in this work are not limited to ferromag- 
netic nanowires on substrates, we have chosen materials, 
geometries and current densities that are typical for spin- 
torque driven domain wall motion studies in ferromag- 
netic nanowires in order to illuminate the role of Joule 
heating in this active research area.^^lti^ 

Secti on [ll| introduces the method underlying the work. 
Section |III| reports results from a number of case stud- 
ies. Starting from the heating of a nanowire of infinite 
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length without constrictions (Sec. Ill A I, we introduce 



a symmetric constriction in a finite-length wire where a 
cuboidal part of material has been removed (Sec. Ill B ) 
to demonstrate the additional heating that results from 
an increased current density in close proximity to the 
constriction and in the constriction (as in Ref . [T2|) . Sec- 
tion Unci stud ies a nanowire with a notch-like constric- 
tion (triangular shape removed from the wire on one 
side only) that is placed on a silicon nitride substrate 
of 100 nm thickness (as in Ref. fT3HT5|) . The same system 
is studied in Sec. IIIIDI where the silicon nitride mem- 
brane is replaced with a silicon wafer with a thickness of 
the order of 500 /im. A zigzag wire on the same silicon 
wafer, as experimentally investigated in Ref. [3l is stud- 
ied in Sec. |III E| We simulate a straight nanowire with- 
out constrictions placed on a diamond substrate as in 
Ref. fT6l in Sec. |IIIFl In Sec. IV A we investigate and dis- 
cuss the applicability of the analytical temperature calcu- 
lation model of Ref. [TTl and deduce an analytical model 
valid for quasi two-dimensional systems such as mem- 



brane substrates in Sec. |IVB[ We briefly discuss free- 
standing and perpendicular nanowires in Sec. |Vj before 
we close with a summary in Sec. |VI[ 



given desired accuracy within their error tolerance set- 
tings. The majority of the other case studies was simu- 
lated using ANSYS. We have taken material parameters 
(see Tab. |l]) appropriate for room temperature, and have 
treated each material parameter as a constant for each 
simulation, i.e. here a temperature dependence of the 
material parameters is not taken into account. 



III. SIMULATION RESULTS 



Case studies 1 and 2 (Sec. Ill A and IIIB) investigate 
a nanowire without a substrate. Case study 3 (Sec. Ill C 
studies a nanowire on (2d) silicon nitride substrate mem- 
brane. Case studies 4 and 5 (Sec. HID and Sec. HIE) in- 



vestigate nanowires on (3d) silicon wafer substrates, and 
case study 6 (Sec. IIIF) reports from a nanowire on a 
(3d) diamond substrate. 

We refer to the substrate as two-dimensional where the 
wire length is much greater than the substrate thickness 
(but still carry out numerical calculations by discretis- 
ing space finely in all three dimensions). We call the 
substrate three-dimensional if the substrate thickness is 
much greater than the wire length. 



II. METHOD 

A current density j(r, t) can be related to the change of 
temperature T(r, t) of a material as a function of position 
r and time t using 



dT 

Ik 



pC 



|.^>V=T + 0) (1) 



with k the thermal conductivity (W/(Km)), p the den- 
sity (kg/m^), C the specific heat capacity (J/(kgK)), 
^ Wx^ + sj/^ + ^ Laplace operator, T the tem- 
perature (K), and Q a heating term (W/m'^). The Joule 
heating of a current density j in an electrical field E is 
given by 
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(2) 



where cr is the electrical conductivity (S/m = l/(f2m)) 
and we have used j = trE. 

Equation ([T]) becomes trivial to solve if we assume a 
uniform current density in a slab of one material with 
constant density, constant thermal conductivity, and con- 
stant specific heat capacity (see Sec. Ill A). In general, 



for samples with geometrical features or spatially inho- 
mogeneous material parameters, the problem becomes 
quite complex and can often only be solved using nu- 
merical methods. For the work presented here we have 
used the simulation software suites ANSYS 12. d^^, Com- 
sol multi-physics,^ and the Nsim multiphysics simula- 
tion libraryiSJ that underpins the Nmag^^ micromagnctic 
simulation package. All three tools were used for case 
study 2 (Sec. IIIB I and produce identical results for a 



A. Case study 1: Uniform current density 

Initially, we study the extreme case of no cooling of 
the ferromagnetic conductor: neither through heat trans- 
fer to the surrounding air, nor to the substrate and the 
contacts. This allows to estimate an upper limit of the 
heating rate and the consequent change in temperature 
over time. 

Assuming a uniform current density j , uniform conduc- 
tivity a and initially uniform temperature distribution T, 
equation ([I]) simplifies to 



dT 
d^ 



.r 

pCn 



(3) 



All the parameters on the right hand side are constant, 
and thus the temperature T will change at a constant 
rate oi p/{pCa). As only the ferromagnetic conducting 
nanowire can store the heat from the Joule heating, the 
temperature has to increase proportionally to the heating 
term Q which is proportional to j^. 

Using material parameters for Permalloy (C ~ 
430 J/(kgK), p = 8700 kg/m^, a = 1/(25 • 10~* f^m) = 
4 • 10^ (r2m)~^, i — 10^^ A/m^), we obtain a change of 
temperature with time 



dT 
d^ 



pCa 



6.683 • 10'" K/s = 66.83 K/ns. (4) 



The parameters for this case study 1 and all other case 
are summarised in Tab. |lj For the permalloy wire, we 
have chosen parameters for case studies 1 to 4 to match 
the experimental work in Ref. [T5HT^ Where possible. 
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constriction (SOnm x 20nm x 20nm), not to scale 
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FIG. 1. Sample geometry for case study 2: slab with constric- 
tion (Sec. |IIIB||. 



parameters measured as part of the experiments have 
been used and have been complemented with literature 
values (see Tab. |l] for details). 

The immediate conclusion from this is that the temper- 
ature of the sample cannot increase by more than 66.83 K 
per nanosecond if the current density of 10^^ A/m^ is not 
exceeded and if the current density is uniform within the 
whole sample for the chosen material parameters. 

A current pulse over 15 nanoseconds has the potential 
to push the temperature up by just over 1000 degrees 
Kelvin, and thus potentially beyond the Curie tempera- 
ture. 

The substrate on which the ferromagnetic conductor 
has been grown will absorb a significant fraction of the 
heat generated in the conductor, and thus reduce the 
effective temperature of the magnetic material. The con- 
tacts play a similar role. On the other hand, any con- 
strictions will result in a locally increased current den- 
sity, which - through the term in Q in Eq. ([2| and 
([T]) - results in significantly increased local heating. We 
study the balance of these additional heating and cooling 
terms in the following sections in detail. 



B. Case study 2: Constrictions 

The effect of a constriction will vary strongly depend- 
ing on the given geometry. The resulting current density 
distribution and spatial and temporal temperature dis- 
tributions are non-trivial. 

Fig. [l] shows the geometry used for this study (as 
in Refs. |32J and 155]) : a bar with dimensions = 
1000 nm, Ly = 50 nm, and Lz = 20 nm. The origin is 
located in the center of the slab, i.e. the two opposite 
corners of the geometry are at [—500, —25, —10] nm and 
[500, 25, 10] nm. The constriction is placed at the center 
of the bar, and reduces the dimensions to ^constrict _ 
20 nm over a length of /v™"'''™' = 50 nm. In compari- 
son to case study 1 (Sec. IIIAl, the current distribution 
is non-uniform in this geometry and therefore the local 
Joule heating and the resulting temperature field will be 
non-uniform. We thus need the thermal conductivity for 
Ni8oFe2o permallojSS] fc = 46.4 W/(K-m) for these calcu- 
lations because the V^T term in equation (nj) is non-zero. 

Fig. [2] shows a temperature profile AT(xj through the 
constricted shape after application of a current density 




FIG. 2. Temperature profile AT{x) in the constricted geom- 
etry shown in Fig. IT] at positions [x,0,0] for t — Ins and a 
current density of 10 ^ A/m^ in the unconstricted ends of the 
slab. 



of lO^^A/m^ at ten different times i as a function of po- 
sition X. We use the notation AT{x) instead of T{x) 
to indicate that this is the change of the temperature T 
relative to the initial temperature, for which we assume 
room temperature (w 300 K). The top thick black line 
shows the temperature distribution after 1 ns, the other 
9 lines show earlier moments in time in successive time 
steps of 0.1 ns. 

The peak around a; = is due to the constriction: the 
reduced cross section (in the y-z plane) results in an in- 
creased current density in the constriction. The Joule 
heating term ([2| — which scales proportional to — 
is increased accordingly in the constriction. During the 
1 ns, diffusion of heat takes place and results in the in- 
crease of temperature outside the constriction, and si- 
multaneously a reduction of the speed of increase of the 
temperature in the constriction. The importance of heat 
diffusion can be seen by the width of the peak around 
X = in Fig. ^ 

The increase of temperature AT after t = 1 ns at the 
end of the nanowire (i.e. x = 500 nm and x — — 500nm) 
is 66.90 K. This is about 0.07 K higher than in the previ- 
ous case study in Sec. [Ill A] where a nanowire without a 
constriction was studied and the corresponding value of 
temperature increase is 66.83 K. This difference of 0.07K 
after 1 ns at the end of the wire originates from the extra 
heating in the constriction around a; = and diffusion of 
this heat through the wire. 

We can also estimate the maximum possible tempera- 
ture increase in the constriction by using equation ^ 
under the assumption that there was no diffusion of 
heat {i.e. fc = 0), and that we obtain a current den- 
sity of j = 2.5 • lO^^A/m^ in the constriction (due to 
the reduced cross section from L.^Ly — 1000 nm^ to 
2^_^^constrict ^ ^^2y rpj^g resulting rate for tem- 
perature increase in the constriction is here 417 K/ns, 
and thus much higher than the maximum temperature 
of 115.6 K that is reached after Ins when taking into ac- 
count the diffusion of heat from the constriction into the 
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Case Study 


Parameter 


Value 


Unit 


Reference 


1-4 


Py 


25 


pf2 cm 


measurecJ^ in experimenlP^'^ R — 280 fl (case study 3) 


5 


Py 


42 


jjbQ, cm 


measured in experimeniP, 7? = 5000 


6 


Py <7-i 


39 


fjbQ, cm 


measured in experimeniP^, R = 675 ft 


1-6 


Py C 


0.43 


J/(gK) 


Ref.Ellp. 252, 6.00cal/(molK) 


1-6 


Pyk 


46.4 


W/(Km) 


Ref.[23]p. 1140, r = 300K 


1-6 


Pyp 


8.7 


g/cw? 


Ref.[HTab. I, lattice constant 3.54 A 


3 


Si3N4 C 


0.7 


J/(gK) 


Ref. EH 


3 


Si3N4 k 


3.2 


W/(Km) 


Ref. |25] 


3 


Si-jNi p 


3 


g/cm^ 


Ref. uni 


4, 5 


Si c 


0.714 


J/(gK) 


Ref.Ei Cp = 20.05 J/(molK) 


4, 5 


Si k 


148 


W/(Km) 


Ref.IWIp. 1-588 


4, 5 


Si p 


2.33 


g/cm^ 


Ref. 121!. Tab. Ill 


6 


Diamond C 


0.53 


J/(gK) 


Ref. EE Cp = 6.37J/(molK) 


6 


Diamond k 


1400 


W/(Km) 


Ref. |30l Fig. 3, Type lb 


6 


Diamond p 


3.51 


g/cm^ 


Ref.EB 



TABLE I. Material parameters used in the simulations (a electric conductivity, a ^ electric resistivity, C specific heat capacity, 
k thermal conductivity, p density). 



unconstricted parts of the nanowire. This comparison 
shows the drastic influence of diffusion of heat on the 
temperatures in the nanowire. 



Case study 3: Nanowire with a notch on a 
silicon nitride substrate membrane 



In this section we investigate a more reahstic system 
following the work by Im et al^^ for which we take into 
account the heat dissipation through the substrate. Im et 
al. study critical external fields for domain wall pinning 
from constrictions, and subsequent workJ^D? study spin- 
torque driven domain wall motion for this geometry. 

Here, we choose the geometry where the domain wall 
de-pinning field was most clearly defined: a permalloy 
nanowire with dimensions = 5000 nm, Ly = 150 nm, 
Lz = 30 nm (see Fig. 1 (b) in Ref. [T3] and lower left 
subplot of Fig. 3 therein) . 

The resistivity of a permalloy thin film strongly de- 
pends on its thickness; in this simulation we scale the 
permalloy resistivity a so that the total resistance of the 
device matches the resis tance value « 280 fl reportecP^' 
for the experiments.'i^^This leads to a resistivity cr"^ = 
25 iiflcm, in line with published data on permalloy resis- 
tivity (Fig. 1(b) in Ref. |M1 t = 30 nm). We use the 



same value for a in case study 1 (Sec. Ill A I, case study 
2 (Sec. Ill B ) and case study 4 (Sec. Ill D ) to allow better 



comparison between the results. 

This nanowire is lithographically defined centrally on 
a silicon nitride membrane that is 100 nm thick (as in 
Ref. I13p . The membrane used in the experiment^'^ was 
purchased from Silson Ltd.'^ According to the man- 
ufacturer the silicon nitride membrane was grown us- 
ing low pressure chemical vapour deposition (LPCVD). 



The thermal properties of LPCVD grown silicon nitride 
films depend on the details of the growth process as can 
be seen i n the range of parameters being cited in the 
literature.l^^'^^1^ We assume values of fc 3.2 W/(mK), 
C = 0.7J/(gK), and p = 3g/cm3 as in Rcf.^S, 

For the modelling of the membrane we use the shape of 
a disk with 0.5 mm radius and 100 nm thickness. Prepa- 
ration of nanostructures on membranes is required for ex- 
periments using synchrotron hght, that has to transmit 
through the sample. Such experiments give access to si- 
multaneous time- and space re solut ion on the nanometer 
and the sub-nanosecond scale .1^^ Perfect thermal con- 
tact between the wire and the silicon nitride substrate 
disk is assumed. The center of the wire contains a 45 nm 
wide triangular notch on one side, and the geometry is 
sketched in Fig. 3(b) A current density of lO^^A/m^ 



is applied at the ends of the wire over a time of 20 ns, 
similar to recent experiments such as in Refs.[Hl ll51 and 

HOI 



Figure p (a) | shows an overview of the geometry and the 
computed temperature distribution after 20 ns. The 30% 
notch is just about visible on the right hand side of the 
nanowire halfway between the ends of the wire. The cut- 
plane shown in the right hand side of Fig. 3(a) as an inset 
shows the temperature distribution in the nanowire in the 
y-z plane in the center of the constriction (as indicated 
by a white line and semi-transparent plane in the main 
plot). Figure 3(b) shows the notch geometry in more 



detail (not to scale). 

Fig. 3(c) shows the temperature profile AT{x) after 
20 ns taken along a line at the top of the nanowire (the 
same data is encoded in the colours of Fig. 3(a) although 
more difficult to read quantitatively). We see that most 
of the wire is at an increased temperature around 270 K. 
The maximum temperature is found at the constriction: 
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FIG. 3. Joule heating in a permalloy nanowire with a notch on a silicon nitride membranj^ as discussed in Sec. 
membrane is modelled as a disk with height 100 nm and radius 0.5mm. The current density is j = 10^^ A/m^; for material 
parameters see Tab.|l] (a) temperature distribution Ar(r) in the nanowire and substrate after t = 20 ns, (b) geometry of the 
model (not to scale) and the plotting path through the nanowire (dashed line) used in Fig. 3(c) (c) temperature profile AT{x) 
after t = 20 ns along the length of the wire (following the plotting path shown in Fig. 3(b) I with a silicon nitride membrane 
substrate, (d) maximum (dotted line) and minimum (dashed line) temperature in the silicon nitride membrane and maximum 
temperature (solid line) in the wire as a function of current pulse length. The solid a nd dot ted line coinci de at this scale. The 
dash-dotted line T^"^ 
further discussion. 



is the prediction of the analytical You model discussed in Sec. IV A See also Sec. IV B and Fig. [t] for 



as in Sec. [Ill 5] the current density is increased here due 
to the reduced cross section in the y-z plane. 

In contrast to the previous example (Sec. IIIB 
and Fig. [2]), the temperature peak at the constriction 
is less pronounced, as a smaller amount of material is 
absent and consequently the current density and the as- 
sociated increase in Joule heating is smaller. 

The maximum temperature increase reaches 290 K. 
The zero level in the simulation corresponds to the tem- 
perature at which the experiment is started: if the ex- 
periment is carried out at 300 K we expect the maximum 
temperature to be 590 K after 20 ns, which is below the 
Curie temperature (~ 840 K) for Permalloy.'^ 



From Figs. 3(a) and 3(c) we can see that the temper- 
ature at the ends of the wire is lower than near the mid- 
dle: for a; < 1 /im and x > 4 fim temperature decreases 



towards « 215 K at the ends of the wire. This is due to 
more efficient cooling through the substrate: at the ends 
of the wire there is substrate to three sides rather than 
two as in the middle parts of the wire. 

The importance of the substrate in coohng the 
nanowire can be seen if we use equation ^ to compute 
the temperature after 20 ns for a wire with the same ge- 
ometry but without the notch and without the substrate: 
the heating rate is dT/dt — p/{pCa) = 66.83K/ns as 
in Sec. Ill A because the geometry does not enter that 
calculation. Without the substrate, we would have a 
temperature increase of ss 1336 K after 20 nano seconds 
(even without taking the extra heating from the notch 
constriction into account). 

In Fig. |3(d)[ the solid line shows the maximum tem- 
perature in the nanowire, the dotted (dashed) line shows 
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the maximum (minimum) value of the temperature in the 
sihcon nitride membrane substrate and the dash-dotted 
line shows data computed using the model by You et al^ 
as a function of time over which the current pulse is ap- 
plied. Note that the You model has not been derived to 
be used for such a thin substrate and that the large de- 
viation for large values of t is thus expected. We discuss 
this in detail in Sec. IIVI 

We can use the maximum wire temperature graph 
to determine the length of the current pulse that can 
be maintained until the temperature is pushed up to 
the Curie temperature, or to the material's evaporation 
temperature (note that the plot shows the temperature 
increase since the start of the experiment, not abso- 
lute temperature). If the experiment is carried out at 
room temperature (w 300 K) , then the Curie tempera- 
ture (w 840 K) is attained with a pulse duration of ap- 
proximately 60 ns. The material starts to melt and evap- 
orate for t > 300 ns, once the melting temperature of 
« 1450 °C is exceeded. 

Asymptotically, the maximum temperature in the wire 
is proportional to the logarithm of pulse duration for a 
2d substrate and a point-like heating source. This regime 
is entered when the nanowire and the constriction have 
been heated up to a steady state, and from there on 
the temperature in the nanowire increases logarithmi- 
cally with time while the heat front propagates from the 
center of the substrate disk towards the disk's boundary. 

The maximum temperature in the substrate (dotted 
line) is assumed at the interface between the nanowire 
and the substrate, in the location where the nanowire is 
hottest. Because of the assumption of perfect thermal 
contact between wire and the membrane substrate, the 
temperature at the top of the membrane is the same as 
the temperature at the bottom of the wire. The difference 
between the maximum wire temperature and maximum 
substrate temperature thus provides an indication for the 
temperature gradient found in the wire. In Fig. |3d| the 
two lines nearly coincide. 

The dashed line in Fig. |3(d)] shows the minimum value 
of the temperature taken across the combined system of 
nanowire and substrate. It starts to deviate from zero 
when the heat front has propagated from the center of 
the silicon nitride substrate disk to the boundary. In 
our example, that happens after w 10^ /iS when the sim- 
ulated temperature increase of the wire would be over 
14,700 K. The asymptotic logarithmic behaviour of the 
temperature increase in the wire is maintained only until 
the heat front reaches the boundary of the substrate. 



than silicon nitride. Furthermore, the extra thickness of 
the wafer gives more space for heat to dissipate. In this 
section, we model the same nanowire with a notch ge- 
ometry as above in case study 3 (Sec. Ill C\ but place 
the nanowire on a silicon half-sphere of radius 0.5 mm 
instead of a 100 nm thin silicon nitride disk membrane. 

Fig. 4(c)| shows the temperature profile of the nanowire 
with notch placed on a silicon wafer after a 20 ns pulse 
(current density at the wire ends is lO^^A/m^). The 
maximum temperature increase is 17 K and should be 
compared with Fig. 3(c) [ where the nanowire was placed 
on a (much thinner) silicon nitride substrate and the 
maximum temperature increase was 290 K. 

Clearly, compared to the (2d) silicon nitride mem- 
brane, the (3d) silicon wafer is much more efficient at 
diffusing heat: the temperature increase is approximately 
a factor 17 smaller. 

Fig. |4(d)| shows the minimum and maximum temper- 
ature for the silicon wafer substrate, and the maximum 



temperature of the nanowire, analog to Fig. 3(d) for the 
nanowire on the silicon nitride membrane. We can iden- 
tify three regimes: (i) for small t the maximum tempera- 
ture in the wire increases approximately proportional to 
the logarithm of time, (ii) For 1 /iS < i < 10, 000 /iS, the 
maximum temperature stays constant in the wire. This 
is the (3d) steady-state regime where the heat front prop- 
agates from the center of the silicon substrate half-sphere 
to the surface of the sphere, (iii) Approximately for 
t > 10, 000 /is, the maximum wire and substrate tempera- 
tures and the minimum substrate temperature start to in- 
crease again simultaneously. This indicates that the heat 
front has reached the surface of the half-sphere shaped 
wafer substrate, and that the heat from the nanowire can- 
not be carried away through the propagating heat front 
anymore. 

The maximum temperature increase of w 19 K is 
reached after ^1/is and is maintained until t — 10 ms, 
when the heat front reaches the boundary of the wafer 
and the whole wafer begins to heat up. The tempera- 
ture gradient from top to bottom of the wire is of the 
order of 5 K (difference of maximum wire and maximum 
substrate temperature) . 

For a system with an infinite 3d substrate we expect 
in general that the nanowire temperature stays constant 
once regime (ii) of the heat front propagation has been 
reached. Regime (iii) would not exist for an infinite 
substrate, or a substrate that is efficiently cooled at its 
boundaries. 



D. Case study 4: Nanowire with a notch on a 
silicon wafer substrate 



A possible way to increase the maximum pulse length is 
to use a silicon wafer instead of the silicon nitride mem- 
brane, which is typically done in high frequency trans- 
port experiments. Silicon is a better conductor of heat 



E. Case study 5: Zigzag-shaped nanowire on a 
silicon wafer substrate 



In this section we investigate a nanowire geometry and 
experimental set up as reported in Ref . |3| where current 
densities of 2.2 x 10^^ A/m^ are applied for 10 /is. 

Figure 5(b)| shows the permalloy nanowire consisting 
of three straight 20 /im segments connected by 45° bends 
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FIG. 4. Joule heating in a permalloy nanowire with a notch on a silicon wafer as discussed in Sec. HID (geometry of wire from 
Refs. [T^and llSf) : (a) temperature distribution Ar(r) in the nanowire and substrate after t = 20 ns, (b) geometry of the model 
(not to scale) and the plotting path (dashed line) used in Fig. 4(c) (c) temperature profile AT{x) after t = 20 ns along the 
length of the wire on a silicon wafer substrate, (d) maximum and minimum temsperatures of the silicon substrate, maximum 
temperature of the wire, and the prediction of the You model T^"^ as a function of pulse length; tc is the characteristic time as 
described in Sec. |IV A| The silicon wafer is modelled as a half sphere (flat side attached to nanowire with notch) with radius 



0.5 mm. The current density is j — 10^^ A/m'^ at the end of the wire. For material parameters see Tab.jl] 



of radius 2 /im (see also Fig. 1 in Ref. |3]). The wire is 
500 nm wide and 10 nm thick, and is placed on a sili- 
con wafer substrate which is modelled as a silicon half- 
sphere of radius 0.5 mm as in the previous example in 
Sec. |III D[ For the simulation we scale the resistivity so 
that the total resistance of the device matches the value 
5 kil reported in Ref. 3 This value is achieved with 
resistivity — 42/iOcm, which agrees well with pub- 
lished data on permalloy resistivity (Fig. 1(b) in Ref. [Ml 
t=10n m). 

Figure [5(a) and 5(c) shows the temperature profile af- 
ter 10 fis. The maximum temperature increase does not 
exceed 133 K despite the large current density and the 
long pulse duration. The ends of the zig-zag wire are sig- 
nificantly cooler (below 80 K) than the center which can 
be attributed to more silicon wafer substrate accessible 
to carry away the heat that emerges from the nanowire. 

Figure 5(d) shows how the maximum wire temperature 



(which is located in the middle segment of the zig-zag 



wire) and substrate minimum and maximum tempera- 
tures change over time. As before, we can identify three 
regimes: (i) initial heating of wire and substrate, (ii) 
steady state with heat front propagating through sub- 
strate and (iii) general heating of wire and substrate 
when heat front reaches the substrate boundary in the 
simulation model. Due to the larger nanowire geometry, 
region (ii) is less pronounced here than in Fig. |4(d)[ 



Case study 6: Straight nanowire on diamond 
substrate 



Recently, realizations of very high current densities 
have been reported when the permalloy nanowire is 
placed on a diamond substrate^ instead of silicon ni- 
tride or silicon. This finding is useful for time-integrating 
experiments where current excit ations exceeding the mi- 
crosecond timescale are requiredj^^^^ In this section, we 
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FIG. 5. Joule heating in zigzag permalloy nanowire on a silicon wafei^^ (Sec. Ill El: (a) temperature distribution Ar(r) in 
the nanowire and substrate after t = 10 /iS, (b) geometry of the model (following Ref. [H]) and the plotting path (dashed) for 
Fig. 5(c) (c) temperature profile AT{x) after t — 10 along the length of the wire (following the plotting path shown in 
Fig. 5(b) I, (d) maximum and minimum temperatures of the silicon substrate, maximum temperature of the wire, and the 
prediction of the You model as a function of pulse length; tc is calculated based on the distance between the opposite ends 



of the wire L — 56 nm. The current density is j 
material parameters see Tab. [T] 



2.2 X 10 A/m ; the resistivity of permalloy a = 42 ^ilcm. For other 



investigate the geometry and experiment described in 
Ref.ie, 



Figure |6(b)| shows a rectangular permalloy nanowire 
with dimensions 25 /zm x 650 nm x 22.5 nm which is grown 
on a diamond crystal substrate. The thermal conductiv- 
ity of diamond depends on the purity of the crystal. For 
the simulation, we conservatively assume a relatively im- 
pure synthetic type lb crystal with thermal conductivity 
k = 1400W/(Km) (Fig. 3 in Ref. W^. A purer crys- 
tal would have higher thermal conductivity and would 
be more efficient at cooling the nanowire. The density 
and heat capacity of diamond are given in Tab. [l] As be- 
fore, perfect thermal contact between the wire and the 
substrate is assumed. In the original experiment ,1^^ the 
whole device was placed in a cryo bath. The area of con- 
tact between the diamond crystal substrate and the bath 
is large compared to the size of the nanowire, therefore 
the temperature difference at the contact layer is likely to 
be small. In the simulation, the effect of the bath can be 



represented by using an infinite medium of diamond. In 
line with the case studies 4 and 5 in the previous sections, 
we use a half-sphere shape for the diamond substrate (ra- 
dius of 0.5 mm). The resistivity of the wire material is 
scaled to 39 /iilcm so that the total resistance of the wire 
is 675 n as in Ref. \M 

Figure 6(a)| and |6(c)| show the temperature profile after 
application of a current density of 1.5 x 10^^ A/m^ for 
1 /iS, and Fig. |6(d)| shows the maximum and minimum 
temperatures in the system as a function of time. 

From Fig. |6(d)] we can see the three regimes as in case 
study 4 and 5: regime (i) shows logarithmic increase 
of temperature with time up to approximately 0.2 fis. 
For larger t the temperature remains constant during 
regime (ii). For t greater than approximately 1000 /is, we 
reach regime (iii) where the heatfront reaches the sub- 
strate boundary. The maximum temperature increase 
in the wire is not exceeding 21 K for times smaller than 
1000 lis, and the substrate temperature increase stays be- 
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FIG. 6. Joule heating in a permalloy nanowire on a diamond crystal substrat^^ (Sec. Ill F I: (a) temperature distribution in the 
nanowire and substrate T after t = 1 /is, (b) geometry of the model and the plotting path, (c) temperature profile Ar(2;) after 
t = 1 us along the plotting path, (d) maximum and minimum temperatures of the diamond substrate, maximum temperature 
of the wire, and the prediction of the You model as a function of pulse length. The current density is j = 1.5 x 10^^ A/m^; the 
resistivity of permalloy a^^ — 39 fiflcm; for other material parameters see Tab. [l] 



low « 16K in that period. The plot shows that the 
precise current density pulse duration is not so critical: 
if t is between 1 fis and 1000 jis the nanowire attains ap- 
proximately the same temperature. It is this steady state 
temperature increase value that should be compared with 
the experiment. 

The difference between the maximum substrate tem- 
perature and the maximum wire temperature reflects 
a temperature gradient from top to bottom in the 
nanowire. This is just about visible in the inset in 



Fig. 6(a 



For larger times than 1000 /iS, we find that the maxi- 
mum and minimum temperature of the substrate and the 
maximum wire temperature start to increase simultane- 
ously. This is an effect of the finite size of the substrate in 
the model, or, equivalently, the lack of the modelling of 
heat transfer away from the diamond substrate through 
the cryo bath. There are other, in comparison to the cryo 
bath less important cooling contributions, such as elec- 
tric contacts, substrate holder, and surrounding gas, that 



have not been considered here. For the interpretation of 
the simulation results for the experiment in Ref. il6\ we 
need to ignore the regime for t > 1000 /is. 



The shape of the temperature profile (Fig. 6(c) ) is very 
similar to the zigzag wire temperature profile (Fig. 5(c) I. 



In the original experiment, a continuous current corre- 
sponding to a current density of 1.5 x 10^^ A/m^ heated 
the wire by about 230 K (Fig. 4 in Ref. [TO)) . 



In our simulation the nanowire changed temperature 
by less than 21 K. A possible explanation for this discrep- 
ancy is that the contact and heat transfer between the 
wire and the diamond was imperfect. 
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IV. ANALYTICAL MODELS 

A. Model by You, Sung, and Joe for a nanowire on 
a (3d) substrate 

You, Sung, and Jo^^^' have provided an analytic ex- 
pression T^'^{t) to approximate the temperature T{t) of 
the current-heated nanowire on a three-dimensional sub- 
strate as a function of time t (Eq. (16) in Ref . [TT|) : 



nka 



(5) 



where w and h are the width and height of the wire, a is 
the wire conductivity, j is the current density, and fc, p, 
and C are the thermal conductivity of the substrate, mass 
density of the substrate, and specific heat capacity of 
the substrate. We have used their adjustable parameter 
a = 0.5 for calculations shown in Figs. [3] to [7] 

We use name T^'^{t) for the equation from You, Sung 
and Joe - which is applicable to 3d-substrates - to em- 
phasize the difference to the simil ar looking equation for 
T^'^(t) that is derived in Sec jlVB and which is applicable 
to 2d-substrates. 

For derivation of Eq. ([s]), it is assumed that the 
nanowire is infinitely long, and attached to a semi-infinite 
substrate. While the thickness h and width w of the wire 
enter the derivation to compute the Joule heating due to 
a given current density, the model does not allow for a 
temperature variation within the nanowire nor does the 
nanowire have a heat capacity in the model. Within this 
model, a heat front of (half-)cylindrical shape (cylinder 
axis aligned with the wire) will propagate within the sub- 
strate when the wire is heated. Thus, there is transla- 
tional invariance along the direction of the wire. 

We start our discussion with the zigzag nanowire as 
shown in Fig. [sj Figure |5(d)| shows the temperature pre- 
diction of the You model as a dash-dotted line. It follows 
the maximum temperature in the substrate very closely 
for times up to approximately 2 fis. 

At short times t below 1 ns we can see the You model 
slightly overestimating the temperature in the nanowire 
in Fig. |5(d)[ As the model does not allow for a finite 
heat capacity of the wire, this is expected. As the heat 
capacity of the nanowire is insignificant in comparison to 
the substrate, this overestimation disappears if sufficient 
heat has been pumped into the system. 

The difference between the maximum temperature in 
the wire and the maximum temperature in the substrate 
comes from a temperature gradient within the wire: the 
maximum temperature in the wire is at the top of the 
wire (which is furthest away from the cooling substrate) 
and the maximum temperature of the substrate is found 
at the top of the substrate just where the wire reaches its 
maximum temperature. Due to the assumption of perfect 
thermal contact, the bottom of the wire is exactly at the 
same temperature as the top of the substrate within the 
model description. 



Since the You model does not allow for a temperature 
gradient within the wire, we expect its temperature pre- 
diction to follow the maximum temperature increase in 
the substrate. This is visible in Fig. 5(d) for i > 2ns. 

Regarding the deviation between the You model and 
the simulation results for t > 2 /is, we need to estab- 
lish whether the required assumptions for the model are 
fulfilled. The You model is derived for an infinitely long 
wire on an infinite substrate, whereas the segments of the 
zigzag wire studied here have finite length. In the initial 
stage of heating, the temperature front in the substrate 
will move away from the wire sections with heat fronts 
aligned parallel with the wire. The heat front forms 
a half-cylinder (for each zigzag segment) whose axis is 
aligned with the wire. As long as the diameter of this 
half cylinder is small relative to the segment length, the 
wire appears locally to be infinitely long and the heat 
front propagates in a direction perpendicular to the wire. 
This is the regime where the You model is applicable, and 
which we have labeled as "regime (i)" in the discussion of 
case studies 4 (Sec. |IIID| to 6 (Sec. |IIIF[ ). When the heat 
front has propagated sufficiently far from the nanowire 
to change its shape from a cylinder surface to a spherical 
surface, the You model is not applicable anymore. This 
happens approximately after t > 2 fis. We have referred 
to the spherical heat front propagation in the discussion 
above as "regime (ii)". 

For the zigzag wire study the agreement of the model 
by You et al.^^ with the simulation is thus very good 
within the time range where the model is applicable. We 
note that the wire is relatively long (20 fim per segment) 
and has no constriction. The You model cannot be ap- 
plied for i > 10 ^s because the finite wire length becomes 
important. 

For the nanowire without a constriction on the dia- 
mond substrate as studied in Sec. IIII Fl and shown in 
Fig. [6j the agreement is similarly good; the You model 
temperature follows the substrate temperature very ac- 
curately up to t « 0.1 fis. For larger t, the model becomes 
inaccurate as the finite length of the wire becomes im- 
portant at that point. 



Figure |4(d)| shows for the nanowire with a notch on 
the silicon substrate that for t > O.l^s the gradients of 
both maximum temperature curves approach zero which 
indicates the onset of regime (ii) and implies that the 
You model cannot be applied for t > 0.1 /iS. For smaller 
t < 0.1 IIS, the You model temperature roughly follows 
the maximum substrate temperature with a maximum 
absolute deviation of less than 3 K. 

We have carried out additional simulations (data not 
shown) which have demonstrated that the maximum sub- 
strate temperature line (dotted line in Fig. |4(d)] ) is shifted 
down by a few degree Kelvin if the notch is removed 
from the geometry. The You model temperature and 
the maximum substrate temperature curves then coin- 
cide for 2 < t < 10 ns. If, furthermore, we increase the 
wire length from 5 fiia to 30 fim, the two curves coincide 
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for 2 < t < 300 ns. 

Both the notch and the relative shortness of the wire 
decrease the accuracy of the prediction of the You model: 
the notch roughly shifts all temperature curves up by a 
few degrees whereas the length of the wire determines 
the time when regime (ii) is entered. 

In contrast to the previous examples the nanowire with 
a notch in Fig. [3] is attached to a relatively thin silicon 
nitride membrane of thickness 100 nni (the silicon and 
diamond substrates for the discussion above are of the 
size order of 500 ^m). The You model should not be 
used in the regime of membranes as the model expects 
an infinite substrate. 

However, one could argue that the You model should 
be applicable for very small t until the heat front emerg- 
ing from the wire has propagated through the 100 nm 
thick substrate. Additional simulations (data not shown) 
reveal that this is the case after t w 0.1 ns. 

We summarize that the You model cannot be expected 
to provide accurate temperature predictions for thin sub- 
strates. The deviation of the You curve in Fig. |3(d) orig- 
inates in the inappropriate application of the model to 
a system with a thin, effectively two-dimensional, sub- 
strate. 

For a 'thick', effectively three-dimensional, substrate 
we find that the applicability of the model is limited by 
the finite length of the wire. The maximum time tc up 
to which the You model is appropriate, can be estimated 
by calculating the characteristic time scale of the heat 
conduction equation ([T]). For a wire of length L, this 
characteristic time tr is 



pC 
k 



(6) 



where fc, p, and C are the thermal conductivity, density 
and specific heat capacity of the substrate material. The 
greater the nanowire length, the longer it takes for the 
heat front to assume spherical shape around the nanowire 
heating source. The greater the heat capacity {pC) and 
the smaller the thermal conductivity, the slower is the 
propagation of the heat front in the substrate. The cur- 
rent density does not enter the equation as it only affects 
the temperature and not the time or length scale. 

Substituting the corresponding parameters for each 
case study, we compute the characteristic time tc for 
the nanowire with a notch in case study 4 (Sec. HID) 
to be 70 ns, for the zigzag nanowire in case study 5 
(Sec. HIE) to be 8.8 /iS, and for the nanowire on dia- 



mond in case study 6 (Sec. IIIFl to be 0.2 /is. We have 
added these values to the figures Fig. |4(d)[|5(d)[ and |6(d)| 
and they are in good agreement with the corresponding 
finite element results. 

The time tc marks the transition from regime (i) to 
regime (ii). 

In summary, we find that the You model provides an 
accurate description of the maximum substrate temper- 



ature if used within its bounds of applicability, i.e. dur- 
ing regime (i) for three-dimensional substrates. In the 
You model, the nanowire has no heat capacity and this 
results in the model slightly overestimating the tempera- 
ture for very small t (visible for example in Fig. |5(d)| for 
t < 1 ns) . The temperature within the nanowire can show 
a gradient (hotter at the top, cooler at the interface to 
the substrate), and the You model computes the smaller 
temperature in the wire. For the studies carried out here 
we find this temperature difference to be less than 10 K in 
all cases although this difference depends on material pa- 
rameters and wire thickness. The You model cannot be 
applied for thin, effectively two-dimensional, substrates 
such as the membrane substrate case study in Sec. |III C| 
and Fig. [H 



B. Analytic expression for nanowire on a 
membrane (2d substrate) 

Equation ^ is valid for pulse durations up to the crit- 
ical duration tc (Eq. ([6])) for effectively three-dimensional 
substrates, i.e. substrates whose thickness is sufficiently 
large so that the heat front in the substrate does not 
reach the substrate boundary for t < tc- This condi- 
tion is fulfilled for the case studies in Sec. |III |III E[ 
and |III F| If the substrate is effectively two-dimensional 
(such as in Sec. IIIC), then Eq. (|5| cannot be applied. 

can be adapted to 



In this section, we show how Eq. 
the 2d case. 

The assumptions made in the derivation of Eq. ([s]) re- 
quire a system of nanowire and substrate that is trans- 
lationally invariant in one direction. In Sec. |III Df[nTF| 
this direction was along the long axis of the wire, and we 
refer to this axis as x, and assume that the height h of 
the wire extends along the z axis. 

For the case of a nanowire on a thin membrane sub- 
strate, we can regard the system as two-dimensional by 
assuming invariance in the perpendicular direction. In 
other words, to apply a modified version of Eq. II), we 
assume that the temperature distribution in the mem- 
brane system is invariant along the z axis. We model the 
nanowire as embedded in the substrate (not grown on top 
of the substrate as in the real system) and imagine an in- 
crease of thickness of both the wire and the substrate 
such that they expand from — oo to -l-oo in z-direction. 
The cooling and heating in each slice (in the x-y plane) 
of the nanowire and substrate system is not affected by 
cooling and heating from slices above and below; it does 
not matter whether we consider only one isolated slice (as 
in the real system) or imagine an infinite stack of slices 
closely packed on top of each other. 

In more detail, we first convert the system of the 
nanowire of height h and substrate of thickness d to a 
2d system of equal height. We increase the height of the 
nanowire by a factor c = d/h so that the wire and the 
substrate are now both of height d (assuming that gener- 
ally d > h but the derivation also holds for d < h). This 
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increases the volume of the wire by a factor of c, and thus 
we wiU have to correct down the heat emerging from the 
wire by the same factor at a later point. 

Second, to obtain translational invariance in the z- 
direction we imagine a stack of such identical 2d systems 
on top of each other. Using the substitutions w — > L and 
h w, we obtain 



1400 - 



T 



2d(t)^^llarcsinh 
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arcsinh 



V 0.5L ) 



(7) 



(8) 



The first fraction in Eq. ([7| is based on Eq. ([s]) and in- 
cludes the substitutions w ^ L and h ^ and we 
also substitute w ^ L in the denominator of the arcsinh 
argument, and use a — 0.5. The second fraction (1/c) 
reduces the temperature increase by c to compensate for 
the increase of heating by the factor c above when we 
increased the thickness of the nanowire to the thickness 
of the substrate. The third fraction (1/2) in Eq. Q is 
a correction because the nanowire is now surrounded by 
substrate in all directions, and not only in one half-space 
as in Eq. ([s]), thus the cooling is twice as efficient and the 
temperature increase is halved. 

Equation Q can be used to compute the maximum 
temperature increase T^'^(i) for a nanowire of length L, 
width w, and height h on a, two-dimensional substrate of 
thickness d. 

In contrast to T^'^(t) there is no upper time limit tc 
for the validity of Eq. (Is]) as the emerging heat- front will 
always stay translationally invariant. 

The comparison of T^'^ {t) with the finite element sim- 
ulation results from case study 3 (Sec. Ill C ) is shown in 
Fig. [7j The overall agreement with the simu ation results 
is good for all times t. The heat capacity of the wire is 
not considered in the model for T^'^(i) which is reflected 
in the overestimation of the temperature in Fig. [7| From 
comparative simulations with different material param- 
eters, we find that the agreement of the two curves is 
better for reduced heat capacity of the wire, and better 
for increased thermal conductivity of the substrate. The 
effect of the finite heat capacity of the wire becomes less 
important for longer current pulses, and the two curves 
in Fig. [b] become closer for larger t (not shown here) . 



V. PERPENDICULAR NANOWIRES 

Recent progress in sample growth has allowed to cre- 
ate free-standing nanowires which are grown perpendic- 
ular to their substrate (for example [44l - l46|) . While it is 
outside the scope of this work to investigate these sys- 
tems in detail, we comment briefly on possible analytical 
approximations. For such a free standing nanowire, the 
analytical expression Q is a good first approximation to 
compute its temperature increase as a function of applied 
current density duration. 
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FIG. 7. Comparison of simulated wire temperature from case 



study 3 (Sec. Ill C I with substrate temperatures obtained us- 
ing the analytical expression T^'^ (t) from Eq. (|8| , as a function 
of pulse length. 



For the temperature increase T^{t) of a perpendicu- 
lar nanowire which is completely embedded in a substrate 
material (such as an AI2O3 matrix), a variation of equa- 
tion ([5| can be employed: 



(9) 



In contrast to the nanowire mounted on top of a half- 
space filling substrate ([s]) as studied in section III the 
substrate is here space-filling, and thus twice as effective 
in cooling the system. 



VI. SUMMARY 

We have carried out detailed numerical simulations of 
the current distribution. Joule heating, and dissipation 
of temperature and heat through the nanowire and the 
substrate for a number of experiments and three different 
substrate types. We find that the silicon nitride mem- 
brane (thickness 100 nm) is the least efficient in cool- 
ing a nanowire that experiences Joule heating. Due to 
the quasi-two-dimensional nature of the membrane, the 
temperature in the nanowire will keep increasing propor- 
tionally to the logarithm of time for longer current pulses 
while the heat front (forming a circle in the membrane 
substrate) propagates away from the nanowire, which is 
located in the center of the heat front circle. 

Using a (effectively three-dimensional) silicon wafer 
substrate instead of the (effectively two-dimensional) sili- 
con nitride membrane, there is a qualitative change: once 
the steady state is entered, the heat front propagates 
in three dimensions and keeps the temperature in the 
heated nanowire constant. In addition to the better cool- 
ing through the three-dimensionality of a silicon wafer, 
the cooling is improved further by the thermal conductiv- 
ity of silicon which is more than one order of magnitude 
greater than that of silicon nitride (see Tab. |l]). 

If we replace silicon in the 3d substrate with diamond, 
the cooling is improved again significantly: diamond's 
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thermal conductivity is about an order of magnitude 
greater than that of silicon (see Tab.|l]). 

In addition to these generic insights, we have worked 
out the temperature increase quantitatively for a series of 
recent experimental publications.^ i3 i6 ^}^q model simu- 
lations presented here show for all of them that the tem- 
perature increase due to the Joule heating did not result 
in the temperature exceeding the Curie temperature. 

We compare these results to the approximating but 
analytical model expression provided by You, Sung, and 
Joe^] and investigate the limits of its applicability. We 
provide an estimate for the characteristic time tc over 
which the You model is valid for three-dimensional sub- 
strates. 

Finally, we provide a new analytical expression that 
allows to compute the temperature for a nanowire on a 
two-dimensional substrate in the presence of Joule heat- 
ing. This expression should be of significant value in 
the design and realization of spin-torque transfer studies 
on membranes where experimenters are often operating 
very close to the Curie temperature or even the melting 



temperature, and where no estimate of the wire's tem- 
perature has been possible so far. 

We provide supplementary online materiaP^ to com- 
pute T^'^(t), and T^'^{t) for other materials and experi- 
mental parameters such as current density and current 
pulse duration. 
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